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LIST  OF  SYMBOLS 

Body  radius,  location  of  the  fin  corner  in 
x-direction  (see  Figure  2). 

Aspect  ratio  of  fin  extended  into  the  body. 

Fin  span,  location  of  the  fin  corner  in 
y-direction  (see  Figure  2) . 

Local  chord.  s 

Reference  chord  | = g j c^(y)dy  j. 

Lift  coefficient  based  on  S. 

Derivative  of  lift  coefficient  lift  curve  slope; 
f3CLl  ^ 

■ bda  - °J based  °"  S' 

Normal  force  coefficient  based  on  S. 

Spanwise  normal  force  distribution  (see  Figure  34). 
Derivative  of  normal  coefficient  norma 1 force  slope; 

■ (^r)a  - o)  based  on  s- 

Surface  pressure  coefficient. 

Surface  pressure  difference  ( = (Cp)lower  - (cp)Upper 

Fractional  location  of  local  panel  in  x-direction 
(see  Figure  2). 

Body  diameter. 

Location  of  local  panel  in  y-direction  (see  Figure  2) . 
Half  width  of  local  fin  panel  (see  Figure  3) . 

Local  chord  at  y . 

Local  panel  chord  (see  Figure  3). 

Number  of  half  body  panels. 

Free  stream  Mach  number. 

Total  number  of  half  fin  panels. 
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Total  number  of  panels  on  fin-body  combination 
( = 2n  + 2M  ) . 

p Aerodynamic  matrix  ( = u/U^  per  unit  y)  . 

Ap/q  Load  distribution  ( = ACp)  . 

q Strength  of  source  per  unit  body  surface  area. 

Q Aerodynamic  matrix  (v/Uro  per  unit  y) . 
r Distance  between  two  points. 

R Aerodynamic  matrix  (w/U^  per  unit  y) • 
s Fin  semi  span. 

S Fin  surface  area  extended  into  body. 

AS  Local  panel  area. 

u Perturbation  velocity  component  in  x-direction. 

Free  stream  velocity. 

v Perturbation  velocity  component  in  y-direction. 

w Perturbation  velocity  component  in  z-direction. 

W Velocity  vector. 

x Body  fixed  x-coordinate  parallel  to  free  stream, 
positive  downstream. 

y Body  fixed  y-coordinate  in  the  spanwise  direction 
of  a fin,  positive  in  right  hand  side  direction 
facing  forward. 

z Body  fixed  z-coordinate , positive  upward. 

ZQ  Fin  location  in  z-direction  mounted  on  the  body, 
a Angle  of  attack. 

8 Prandtl-Glauert  parameter  ( = Vl-Mw  ) . 

y Circulation  strength  of  fin  vortex. 
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Circulation  of  fin  vortex. 

Cant  angle  of  fin,  deflection  angle  of  flap  or 
elevon. 

Azimuth  angle  of  body,  positive  counterclockwise 
facing  forward  (see  Figure  4) 

Azimuth  angle  of  body,  positive  clockwise  facing 
forward  (see  Figure  13b) . 

Gradient  of  bound  vortex,  (see  Figure  3),  taper 
ratio. 

Sweptback  angle  at  specified  chordwise  position. 
Apparent  angle  of  attack. 

Distance  between  two  points. 


Body 

Due  to  bound  vortex  (see  Appendix  1) . 
Control  point  on  panel  (see  Figure  3) 
Elevon. 

Fin. 

Flap . 


Local  panel  name 


Normal  to  body  surface. 

Due  to  port  side  vortex  (see  Appendix  1) . 

Due  to  starboard  side  vortex  (see  Appendix  1) 
Tangential  to  body  surface. 

Bound  vortex  (see  Figure  3) . 
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Aerodynamic  performance  of  various  kinds  of  fin-body  combi- 
nations in  inviscid,  subsonic  flow  at  small  angles  of  attack  can 
be  analyzed  by  the  "singularity  method"  (Refs.  1 and  5).  Koerner's 
approach  (Ref.  1)  is  followed  in  this  analysis.  In  general,  the 
singularity  method  has  an  advantage  over  the  slender  body  theory 
in  the  sense  of  less  restriction  in  the  fin  shape.  By  means  of  a 
modern  high  speed  electronic  computer,  the  pressure  coefficient 
can  be  calculated  rather  simply. 

The  surface  pressure  distribution  on  the  rectangular,  the 
sweptback,  and  the  tapered  sweptback  fins  and  fin-body  combinations 
have  been  analyzed.  The  effect  of  compressibility  is  taken  into 
account  by  Goethert's  rule  which  showed  very  satisfactory  with  the 
experimental  data  up  to  a high  subsonic  free  stream  Mach  number. 

For  all  of  these  no  shock  waves  oil  fins  and  fin-body  combinations 
are  treated. 

The  present  analysis  has  been  extended  to  include  the  slender, 
cruciform,  canted  delta  fin  and  fin-body  combination  and  also 
includes  the  case  of  slender,  cruciform,  clipped  delta  fin-bodv 
combination  geometry  v/ith  a deflected  small  control  fin.  The 
results  nave  been  compared  with  the  slender  body  theory  (Ref.  2) 
and  other  singularity  methods  (Refs.  5 and  8). 

The  preset t method,  with  slight  modification,  can  be  applied 
also  to  a fin-body  combination  geometry  of  a fin  v/ith  twist,  camber, 
thickness,  and  arbitrary  plan  form.  From  the  analytical  study, 
it  is  also  clear  that  a viscous  flow  study  is  needed  in  estimating 
the  more  accurate  pressure  distribution  in  the  region  of  the  fin- 
body  junction  as  well  as  the  fin  tip,  even  at  small  angles  of 
attack. 

For  more  condensed  contents  of  the  present  report  and  its 
interaction  to  the  exhaust  jet  plume,  see  Ref.  (12). 
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SECTION  II.  BASIC  ANALYSIS 

In  estimating  the  aerodynamic  performance  on  a fin-body 
combination  geometry,  the  method  based  on  slender  body  theory 
has  been  developed  and  easily  applied.  Very  reasonable  results 
have  been  obtained  (see  Ref.  2).  The  practical  difficulty  in 
applying  the  slender  body  theory  is  in  its  limitation  to  the 
very  slender  fin  configuration  (e.g.,  one  needs  to  assume  that 
/R  <<  1 theoretically,  or  /R  < 1 even  for  a practical  application  to 
be  held.  On  the  other  hand,  the  Prandtl's  well-known  "lifting 
line  theory"  has  approximated  very  well  the  lift  distribution  for 
a rectangular  fin  with  very  large  aspect  ratio.  The  application 
to  the  other  fin  configurations  (such  as  a finite  fin)  seems 
not  that  good  as  far  as  the  accuracy  is  concerned  for  the  surface 
pressure  distribution.  To  overcome  this  deficiency,  Lawrence 
(Ref.  A ) has  developed  the  approximate  solution  of  the  lifting 
line  theory  for  low  aspect  ratio  wings (including  the  rectangular 
shape)  and  delta  fins  with  the  aspect  ratio  from  zero  to  four. 

However,  the  application  of  this  technique  to  the  fin-body  combi- 
nations with  the  other  fin  configurations  seems  to  be  still  in 
difficulty.  A "singularity  method"  is  a method  of  numerical  anal- 
ysis. From  the  study  of  literature,  it  can  be  said  that  such  a 
numerical  approach  will  provide  the  best  estimation  of  the  pressure 
distribution  in  nearly  all  the  fin-body  configurations  with  a good 
accuracy . 

1.  Basic  Consideration 

For  simplicity,  the  following  assumptions  were  made  in 
this  analysis,  i.e., 

(i)  A steady,  inviscid,  uniform,  subsonic  free  stream. 

(ii)  A flat  fin  with  a straight  leading  and  trailing 
edge,  and  a streamwise  fin  tip. 

(iii)  Fins  are  attached  to  a circular  cylindrical  body. 

(iv)  A small  angle  of  attack. 

Based  on  the  singularity  method,  constant  strength  vortex 
singularities  are  used  to  replace  a lifting  fin  on  a horizontal 
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plane.  The  image  vortices  are  placed  inside  the  body  with 
respect  to  the  body  surface  in  order  to  compensate  the  body 
displacement  effect.  To  obtain  a good  aerodynamic  interaction 
for  a fin-body  geometry,  the  source-sink  singularities  are 
distributed  on  a body  surface  and  the  doublet  singularities  are 
distributed  along  the  body  axis.  The  doublet  distribution  gives 
the  effect  of  a body  inclination  to  a free  stream.  The  iteration 
scheme  is  used  so  that  all  the  simulated  lifting  surfaces  can 
eventually  satisfy  the  flow  boundary  conditions  everywhere  on  the 
surface  of  the  fin-body  combination. 

Figure  1 shows  a schematic  model  of  a basic  fin-body 
combination  for  this  analysis.  A fin  (or  wing)  has  been  divided 
to  many  panels  and  each  one  has  been  replaced  by  a horsesWbe -type 
vortex  which  consists  of  three  vortices  (i.e.,  one  bound  vortex 
and  two  free  vortices) . The  body  has  been  divided  also  into  many 
panels,  each  of  which  consists  of  a correct  part  of  a body,  and 
it  has  been  replaced  by  a source  or  a sink  at  its  surface.  (Also, 
set  Figure  51  in  Appendix  1) . 


2.  Basic  Equations  and  Coordinates 

Figure  2 shows  the  coordinate  of  a fin  part,  which 
is  divided  by  a number  of  trapezoidal  panels  which  have  the  sides 
parallel  to  a free  stream.  Each  panel  has  three  vortices  with  a 
constant  strength,  i.e.,  one  bound  vortex,  which  is  located  at 
one  quarter  local  panel  chord  length,  and  two  free  vortices  which 
are  located  at  local  panel  edges,  parallel  to  the  free  stream. 
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The  location  of  a control  point  at  which  the  boundary  condition 
will  be  satisfied  is  assumed  at  three-quarters  of  a local  panel 
chord  length  in  the  center  line  of  a local  panel  as  shown  in 
Figure  3.  In  addition  to  this,  the  Goethert  compressible 
similarity  rule  has  been  employed  to  the  present  analysis,  i.e., 
the  coordinates  normal  to  a free  stream  have  been  multiplied  by 
Prandtl-Glauert  Parameter  6. 

The  following  five  steps  are  essential  to  the  analysis. 

The  steps  were  established  by  Ref.  1. 

(i)  Determining  circulation  distribution  (y^)  on  a fin. 

Normalized  perturbation  velocity  components 
at  any  control  point  on  panel  i induced  by  any  panel  j which  makes 
four  discrete  horseshoe  vortices  (two  are  from  a pair  of  fins 
themselves,  and  two  are  from  their  body  images)  can  be  written  as. 
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where:  u^,  v^  w^  are  perturbation  velocity  components  and  , 

and  are  the  aerodynamic  matrices  given  in  Appendix  1. 

The  third  equation  as  stated  above  is  used  for 
satisfying  the  boundary  conditions  of  any  panel,  i.e., 


(ii)  Determining  Source  Distribution  (qy)  on  the  Body. 

Induced  velocity  components  at  any  point  on 
panel  v on  the  body  surface  due  to  the  fin  and  its  image  vortices 
on  the  control  point  of  panel  j is  expressible  as, 


V 

Sbv 


v = 1,  2M 

where  Fb,  Qb  and  Rb  are  the  same  aerodynamic  matrices  as  those 
in  Equation  (1),  except  taking  the  control  points  on  the  body,  which 
are  located  at  the  center  of  local  panels  (see  Figure  51  in  Appendix 
1) , in  the  former  case. 

In  cylindrical  coordinates , 

vn^  = vb  cos  9 + wb  sin  9 


vt  = -vb  sin  0 + wb  cos  0 

Geometrical  relationships  are  illustrated  in  Figures  4 and  5. 

The  resultant  velocity  component  normal  to  the  body  surface  due  to 
the  fin  vortices  and  the  sources  must  be  varnished,  i.e., 

vnl  + vn2  = 0 <5> 

where  vn2  is  the  velocity  component  normal  to  the  body  surface 
due  to  the  sources  distributed  on  it,  vn2  can  be  written  as, 

v . f -afr'.O.  n cm  (o-q)  a{  d0,  dx,  (6) 

7 i L-n  n*  S 


in  incompressible  flow,  where 


(x  - x')2  + (y  - y')2  + (z  ~ z’)2 


A 


The  integral  equation  (6)  can  be  solved  by  the  iteration  scheme; 
For  the  first  approximation, 


or. 


<ll  (x,0) 


q2  (x,6) 

00  r2*  vnl  (x-  ,e ' ) 

o " 


-2  vnl  (x.e) 

-2  vnl  (x,0)  + 

{1  - cos  (0-9 ' ) } a^ 


(8) 


de'dx'  (9) 
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The  nr  approximation  of  the  source  strength  at  a point  v 

on  a body  surface  which  cancels  the  induced  velocity  v , in 

n J. 

Equation  (4)  normal  to  a body  surface  can  be  written  as 

q - -2  (V  ,)  - 2ZM 

™ ' nl,m“  v-1  2 , 


{ 1-cos  (ep- ev) > 0 a A sv 
TCVV2  + e2  { 1_cos  (Vev)}  fl2] 


Sv  = 0 a A0  Axv 

(iii)  Determining  Downwash  (w^  ) on  Fin  Due  to 
Sources  Distributed  on  the  Body  Surface 


The  downwash  on  panel  u due  to  source  on  panel 


v can  be  written  as, 

W 2M  % 

Wbu  = E t- 


0(zo-zv)  ASV 


K [ <vxv>2+6  { <vyv)2+(zo-zv>2}  ] 3/ 


(iv)  Determining  Circulation  Distribution  (y^) 
for  Second  Iteration 

By  the  downwash  obtained  in  Equation  (12), 
the  corresponding  angle  of  attack  is  induced  on  a fin,  which  can 
be  written  as  W^/U^.  Then,  this  term  should  be  added  to  the 
original,  i.e.,  geometrical  angle  of  attack  in  Equation  (2),  thus 
must  be  written  as, 
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After  has  been  obtained  on  fin  panels, 
the  same  procedure  should  be  repeated  until  y±  will  not  change 
greatly. 

(v)  In  Case  of  Body  with  Angle  of  Attack. 

Downwash  on  a fin  due  to  the  doublet  along 
a body  axis  (as  shown  in  Fig.  6)  can  be  written  as 


2 , 2 y2< 

IU  “ba  (y2+z2)2 


(14) 


This  downwash  has  been  assumed  constant  at  any  spanwise  location, 
and  it  changes  an  apparent  angle  of  attack  (wi)  as  mentioned 
above.  Thus,  in  this  case,  should  be  replaced  by, 


-u».  = 


W,  . W-. 

-ac.  + + _li 

Uoo  U«, 


(15) 


The  aerodynamic  performance  coefficient  can  be  computed  in  the 
usual  manner  for  the  pressure  coefficient,  Cp , on  panel; 


Cpj  . 


— for  fin, 


2u 


U 


lhl 


CO 


for  body. 


> (16) 


where  "+"  is  for  an  upper  surface  of  a fin  and  is  for  a lower 
surface. 
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For  the  load  distribution,  AC  , 

P 


6 <*drv!> 


for  fin, 


(Cpj> lower  - <Cpj> 
surface 


PJ ^ upper  •••  for  body- 
surface 


For  the  spanwise  lift  distribution,  CL(y) , 


n / \ 2b  " 

L y 02l(y)  i=l  Yi 


where  n'  is  a total  number  of  panels  along  the  chord  l(y)  at 
a fixed  y station. 

For  the  lift  coefficient,  CL> 


£ C . AS- 


cT  = izL_ 

L p 


PJ  J 


E AS; 


where  AS..  = 40hj  (x^-x  .)  for  the  fin. 


= 0aA0AXjSin  0j  for  the  body. 


(projected  area  on  horizontal  plane) 


SECTION  III.  RESULTS  AND  DISCUSSIONS  ON  VARIOUS  FINS. 

1.  Rectangular  and  Constant  Chord  Sweptback  Fins. 

As  the  first  example  by  the  singularity  method,  the 
spanwise  lift  distribution  at  Mra  = 0 was  computed  for  the  fin 
with  the  three  different  leading  edge  sweptback  angles,  that  is, 
zero  (rectangular  fin),  30  and  45  degrees.  The  aspect  ratio 
of  these  fins  is  taken  as  all  the  same  equals  to  six.  The  fin 
configuration  and  the  computed  results  are  shown  in  Fig.  7.  This 
figure  shows  the  spanwise  lift  distribution  per  unit  angle  of 
attack  over  the  semi  span  of  the  fin.  The  half  of  the  fin  has 
been  divided  into  the  sixteen  trapezoidal  panels,  i.e.,  four 
equidistantly  in  the  chordwise  direction,  and  four  with  cosine 
in  the  spanwise  direction.  The  present  calculations  have  been 
compared  with  those  in  Ref.  (1),  on  which  the  present  analysis 
was  based.  The  results  are  shown  in  Fig.  7.  It  can  be  seen  that 
there  is  little  difference  between  them.  As  the  ordinate  in 
Fig.  7 is  considered  proportional  to  the  circulation  distribution 
per  unit  angle  of  attack  over  the  semi-span  of  the  fin  (see 
curve  A),  i.e.,  the  rectangular  fin,  is  very  close  to  the  ellip- 
tics.  With  increasing  the  sweptback  angles  of  a fin,  the  lift 
distribution  around  the  center  of  the  fin  is  decreasing,  or  the 
maximum  point  of  the  lift  distribution  moves  to  the  fin  tip,  and 
thus  the  lift  of  the  fin  decreases  accordingly.  This  means  that 
the  characteristics  of  the  fin  will  be  close  to  the  one  of  the 
"delta"  fin  (compare  it  with  Figs.  11  and  12). 

The  typical  chordwise  pressure  distribution  on  a rectangular 
fin  at  the  different  spanwise  positions  is  shown  in  Fig.  8.  The 
same  rectangular  fin  configuration  as  the  one  in  Fig.  7 was  used 
here.  The  pressure  goes  to  infinity  at  the  leading  edge  of  the 
fin  because  of  the  "leading  edge  singularity"  of  the  fin  with 
zero  thickness.  It  can  be  also  been  that  the  Kutta's  conditions 
at  subsonic  speed  are  certainly  satisfied  at  the  trailing  edge  of 
the  fin. 


2.  Compressibility  Effect  on  Pressure  Coefficient  of 
Rectangular  Fin. 

The  variation  of  pressure  coefficient  with  the  free 
stream  Mach  number  at  one  position  on  the  rectangular  fin  is 
shown  in  Fig.  9.  Goethert's  rule  was  used  for  compensating 
the  compressible  effect  in  the  present  analysis  as  mentioned 
before.  The  computed  results  by  using  the  other  compressible 
similarity  rules,  such  as  Prandtl-Glauert  and  KArmdn-Tsien  rules, 
are  also  included  in  this  figure  for  comparisons.  The  Cps 
corrected  by  using  Prandtl-Glauert  and  Karman-Tsien  rules  seem 
to  be  overes timt ted.  The  correction  of  compressibility  based  on 
Goethert's  rule  predicts  satisfactorily  well  for  a fin  with  low- 
aspect  ratio  up  to  a relatively  high  Mach  number  (for  example, 
see  Chapter  13  in  Ref . (10)  ).  With  increase  in  the  aspect  ratio 
of  a fin,  the  difference  of  the  lift  coefficients  computed  by 
using  the  Prandtl-Glauert  and  the  Goethert  rules  becomes  small. 
Fig,  10  shows  the  calculated  example  on  the  lift  coefficient,  for 
a rectangular  fin  with  aspect  ratio  of  twenty.  It  should  be 
noted  that  both  of  the  calculated  results  have  little  difference 
up  to  near  a Mach  number  of  0.8. 

3.  Clipped  Delta  Fin. 

The  pressure  distribution  for  a clipped  delta  fin  con- 
figuration has  been  computed.  The  result  is  shown  in  Fig.  11. 

The  solid  lines  indicate  the  chordwise  pressure  distribution  along 
the  corresponding  chord  lines,  and  the  dotted  lines  indicate  the 
spanwise  pressure  distribution  along  the  one-quarter  chordwise 
length  of  the  local  panel.  As  can  be  seen,  the  clipped  fin  tip 
effect  shows  up  after  the  entire  mid  chord  length  or  toward  the 
trailing  edge.  It  is  relatively  not  affected  in  the  front  portion 
of  the  chord.  The  free  stream  Mach  number  for  this  computed  case 
is  0.7,  the  leading  edge  sweptback  angle  is  77.2°,  and  the  angle 
of  attack  is  taken  as  one  degree.  From  this  figure,  it  can  be 
seen  that  the.  aerodynamic  characteristics  of  the  clipped  delta 
fin  has  the  one  combined  with  those  of  the  sweptback  and 
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Figure  9.  Cp  at  one  point  on  rectangular  fin  changing  with  free  stream 
Mach  number  computed  by  different  subsonic  compressible  similarity 
rules. 
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Figure  10.  Comparison  of  lift  coefficient  of  rectangular  fin  with  large 
aspect  ratio  by  two  different  similarity  rules. 
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Figure  11.  Pressure  difference 
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Fig.  12  shows  the  spanwise  lift  distribution,  and 
that  of  the  delta  fin,  the  geometry  of  which  is  shown  in  Fig.  11, 
is  also  included  in  this  figure  for  a comparison.  Although  the 
actual  positions  in  the  spanwise  direction  are  slightly  different 
in  both  cases,  the  decrease  of  the  spanwise  lift  distribution 
caused  by  a clipping  of  the  fin  tip  can  be  seen  very  clearly  near 
the  fin  tip. 
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SECTION  IV. 


RESULTS  AND  DISCUSSIONS  ON 
VARIOUS  FIN-BODY  COMBINATIONS. 


1.  Rectangular  Fin-Body  Combination. 

The  spanwise  lift  distributions  have  been  calculated 
for  a rectangular  fin  alone  and  a fin-body  combination 
configuration  at  Mm  = 0.  They  agree  well  with  the  original 
calculation  of  Koerner  (Ref.  1)  as  shown  in  Fig.  13a.  The 
discussion  on  the  compressibility  effect,  will  be  given  later. 

The  body  treated  has  a zero  angle  of  incidence,  and  a fin  is 
canted  by  one  degree  (clearly,  for  a more  cant  fin,  a simple 
linear  multiplication  of  amplitude  will  yield  a result).  The 
presence  of  the  body  was  simulated  by  image  vortices  of  the  fin 
and  source-sink  distributed  on  the  body  surface  as  discussed. 

The  computed  results  for  the  fin-body  interaction  agreed  well 
with  a limited  experimental  data  as  shown  in  Fig.  13a,  except  in 
the  region  very  close  to  the  fin  root.  Koerner  pointed  out  that 
this  disagreement  resulted  from  the  presence  of  boundary  layer 
in  a real  flow  on  the  fin  and  the  body  (Ref.  1).  This  suggests 
that  the  viscous  study  is  necessary  in  order  to  obtain  a more 
accurate  result  in  this  small  region.  As  it  can  be  seen  from 
Fig.  13a  , the  mutual  interaction  between  the  body  and  the 
fin  is  rather  strong  (comparing  cases  with  and  without  the  body 
presence) . 

The  chordwise  load  distributions  on  the  body  in  the 
same  fin-body  geometry  as  shown  in  Fig.  13a  were  computed  at 
three  different  azimuthal  positions.  The  results  are  shown  in 
Fig.  13b.  It  can  be  seen  that  the  influence  of  the  fin  on  the 
body  becomes  stronger  as  closer  to  the  fin  root.  It  affects  both 
the  upstream  and  downstream  flow  beyond  the  fin  location.  The 
calculated  chordwise  lift  distribution  on  the  fin  very  close  to 
the  fin  root  is  also  included  in  Fig.  13b  for  a comparison. 

Fig.  14  shows  the  chordwise  load  distribution  on  the 
rectangular  fin-body  combination,  in  which  the  fin  is  not  canted, 


24 


V 


Figure  14.  Longitudinal  pressure  distribution  on  rectangular  fin-body 
combination  (a,  = 6.8#,  M =0). 

D 00 


Figure  16b.  Longitudinal  pressure  distribution  on  body. 


3.  Tapered,  Sweptback  Fin-Body  Combination. 

The  sweptback  fin-body  combination,  with  two  different 
taper  ratio  fins  as  used  by  Koerner  (Ref.  1)  has  been  investi- 
gated. Fig.  17  shows  the  computed  spanwise  lift  distribution  on 
such  a geometry  with  a taper  ratio  X = 1/1.  (where  1-  is  the 
maximum  chord  length  through  the  body,  and  1 is  a fin  tip  chord 
length)  of  zero  and  one  third.  The  sweptback  angle  for  both 
cases  was  taken  30°  at  one  quarter  of  the  local  chord  length.  The 
angle  of  incidence  of  the  body  was  assumed  one  degree  (again,  for 
a larger  angle  case,  a linear  multiplication  of  amplitude  is 
sufficient)  and  the  fin  was  assumed  with  a zero  cant  angle.  The 
lift  distribution  for  X = 0 agrees  well  with  the  characteristics 
of  the  delta  fin-body  combination,  in  spite  of  the  fact  that  it 
is  not  quite  a delta  fin.  Also,  compare  it  with  the  result  in 
Fig.  19a,  a case  of  slender  delta  fin-body  is  shown  although  the 
fin  configuration  is  not  the  same. 

The  spanwise  lift  distribution  of  the  trapezoidal  fin-body 
combination,  the  configuration  of  which  is  taken  from  Ref.  (7) , 
is  shown  in  Fig.  18.  The  trend  of  the  lift  distribution  curve 
has  not  much  difference  with  the  one  in  the  case  of  X = 1/3  as 
mentioned  before,  or  the  fin  part  of  the  clipped  delta  fin  case 
as  shown  in  Fig.  12. 

4.  Cruciform  Slender  Delta  Fin-Body  Combination. 

Fig.  19a  shows  the  spanwise  lift  distribution  on  three 
different  cruciform  fin-body  combinations.  In  the  case  of  the 
cruciform  cant  fin-body  combination,  a pair  of  horizontal  fins  is 
canted  by  three  degrees  so  that  it  can  rotate  about  the  body  axis 
(i.e.,  x-axis)  to  a negative  direction  in  the  sense  of  a commonly- 
accepted  sign  convention.  The  case  of  the  fin  only  yielded  the 
largest  lift  distribution  as  expected.  The  effect  of  the  body 
can  be  seen  clearly  in  this  case  also.  The  cant  fin  case  showed 
the  lowest  lift  distribution  as  expected  (about  this,  for  example, 
see  Ref.  (2)  ) . 


Figure  19b.  Longitudinal  pressure  distribution  on  body  with  cruciform 
slender  delta  fin  (horizontal  fins  are  canted,  6 = -3°). 


Fig.  19b  shows  the  calculated  longitudinal 
pressure  distribution  on  the  body.  It  should  be  noted  that  the 
rate  of  perturbation  pressure  change  is  very  sharp  at  the  leading 
and  the  trailing  edges  with  lower  magnitude  pressure  as  compared 
with  the  large  aspect  ratio  fin  case  (see  Fig.  13b).  This  is 
expected  because  the  fin  is  slender  so  the  influence  of  the  fin 
on  the  body  should  be  small.  This  result  agrees  well  with  the 
concept  of  the  slender-body  theory  which  states  that  the  aerodynamic 
disturbance  influences  only  in  the  cross  plane  to  the  free  stream 
(see  Ref.  (2)  ) . 

5.  Discussions  on  Compressibility  Effect  on  Fin-Body  Combination. 

The  effect  of  the  compressibility  on  a lift 
coefficient  of  a cruciform  slender  delta  fin-body  combination 
is  shown  in  Fig.  20.  (Only  the  angle  of  attack  influence  is  con- 
sidered here  for  the  sake  of  discussion.)  The  fin  area  extended 
through  the  body  as  shown  on  a dotted  line  was  used  as  a reference 
area  for  the  lift  coefficient  computation.  In  general,  as  mentioned 
before,  the  Prandtl-Glauert  rule  resulted  in  an  over  correction 
for  such  a small  aspect  ratio  fin  configuration  (see  Fig.  20). 

It  is  more  appropriate  to  use  the  Goethert's  rule 
to  calculate  the  compressibility  effect  in  a fin-body  combination 
geometry.  For  a slender  fin-body  combination  case,  after  using 
the  Goethert  rule  correction,  the  result  (see  Fig.  20)  agreed  very 
well  with  the  result  computed  by  the  slender  body  theory.  In  view 
of  a wider  range  of  Mach  number  applicability  of  the  slender  body 
theory,  such  an  agreement  is  not  totally  surprising.  In  the  same 
figure,  a result  by  Nielsen's  estimation  (Ref.  11)  is  also 
included  for  a comparison.  For  a non-slender  fin  case,  the 
result  by  using  Goethert's  rule  correction  is  shown  in  Figs.  21 
and  22. 

Fig.  23  shows  the  lift  coefficient  varying  with 
the  freestream  Mach  number  on  the  non-slender  rectangu. ar  fin- 
body  combination.  The  configuration  is  given  in  Fig.  22.  As  it 
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Figure  21.  Local  lift  coefficient  on  fin  body  junction 
showing  effect  of  Mach  number  change,  -30°  swept  fin 
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can  be  seen  very  clearly  from  Figs.  20  and  23,  the  compressi- 
bility effect  of  lift  coefficient  on  the  slender  fin-body 
combination  is  much  more  "blunt",  or  little  change  up  to  quite 
high  subsonic  speed  than  the  one  of  the  non-slender  fin-body 
combination  case. 

6.  Effect  of  Small  Control  Fin  on  Slender  Fin-Body  Combination. 

One  of  the  advantages  of  the  singularity  method  is 
that  the  effect  of  a control  fin  can  be  computed  very  easily. 

The  chordwise  pressure  distribution  has  been  computed  for  the 
slender  delta  fin-  (one  part  of  which  is  used  as  the  control  fin) 
body  combination  at  = 0,  and  is  given  in  Fig.  24.  The  hinge 
line  of  the  control  fin  is  assumed  to  be  located  at  the  three 
quarter  local  chord  length  of  the  delta  fin.  Angle  of  deflection 
of  this  fin  is  three  degrees  in  the  sense  of  the  trailing  edge 
of  the  control  fin  downwards.  The  main  fin  is  assumed  to  position 
at  zero  angle  of  attack.  Although  such  a computational  configura- 
tion of  the  control  fin  position  has  not  a great  meaning  in  a 
practical  application,  yet  the  effect  of  its  controllability  can 
be  seen  very  clearly.  The  hinge  line  can  be  considered  as  a 
singular  point.  And  thus  the  local  pressure  goes  to  infinity. 
However,  this  is  not  a realistic  case.  The  pressure  at  this  point 
cannot  be  infinity  because  of  the  viscous  effect  in  a real  situation. 
It  can  be  seen  that  the  deflection  of  a control  fin  influences 
ever  the  whole  surface  of  the  main  fin.  Such  an  influence  is  very 
strong  near  the  fin  tip.  (Compare  the  present  calculation  with 
that  in  Ref.  9,  for  example.) 
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SECTION  V.  KOERNER’S  APPROACH  COMPARED  WITH 
OTHER  METHODS  AND  DATA. 


1.  Comparison  with  the  Slender-Body  Theory. 

A comparison  of  computed  results  by  the  singularity 
method  and  the  slender-body  theory  had  been  reported  in  the  pre- 
vious subsection.  A surprisingly  good  agreement  by  both  methods 
was  obtained  for  the  lift  coefficient  of  a cruciform  slender 
delta  fin-body  combination. 

The  pressure  and  the  load  distributions  on  a fin-body 
combination  of  the  same  computational  model  as  used  in  the  slender 
body  theory  example  (see  Ref.  (2))  was  calculated.  The  results 
are  shown  in  Figs.  25  and  26.  Fig.  25a  shows  the  pressure  distri- 
bution on  one  horizontal  fin.  It  is  worthwhile  to  note  that  the 
linearized  Cp  expression  was  used  in  the  present  analysis  while 
the  Cp  expression  in  slender-body  theory  was  considered  up  to 
the  second  order  terms.  The  agreement  cf  both  methods  is  excellent. 
The  pressure  coefficient  was  computed  azimuthally , at  the  six  posi- 
tions on  the  half  body,  and  is  shown  in  Fig.  25b.  Note  that  the 

difference  in  the  C expression  is  more  sensitive  on  the  fin  than 

P 

on  the  body.  This  means  that  the  second  order  expression  of  the 
Cp  is  more  important  for  such  slender  fin-body  combination  geom- 
etry even  though  its  absolute  value  is  small.  Fig.  26  shows  the 
load  distribution  on  a horizontal  and  a vertical  fin.  It  can  be 
seen  from  Figs.  25  and  26  that  the  results  of  the  present  analysis 
and  slender-body  theory  have  agreed  satisfactorily  with  each 
other.  From  these  comparisons,  it  can  be  concluded  tl  t both 
methods  will  yield  a similar  result  for  a small  aspect  ratio  fin- 
body  combination  geometry.  This  merely  proves  that  the  singularity 
method  can  be  applied  for  fins  with  a wide  range  of  aspect  ratios 
and,  therefore,  is  more  versatile. 
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2.  Comparisons  with  FLEXSTAB  Singularity  Distribution 
Method  and  Data. 

FLEXSTAB' s method  (Ref.  8)  has  been  applied  to  a 
wing-body  combination  of  an  aircraft  flying  at  subsonic,  tran- 
sonic and  supersonic  speeds  (Refs.  5 and  8).  A comparison  of 
computations  based  on  the  Koerner ' s scheme,  as  adapted  in  this 
analysis,  was  made  to  those  configurations  calculated  based  on 
FLEXSTAB’ s technique.  The  computational  wing-body  configuration 
which  is  taken  from  Ref.  (8)  is  shown  in  Fig.  27,  i.e.,  the  aspect 
ratio  of  the  wing  is  1.65,  taper  ratio  is  0.10,  and  sweptback 
angle  of  the  leading  edge  is  71.2°,  respectively. 

The  paneling  scheme  employed  for  the  present  calcula- 
tions as  well  as  those  by  the  FLEXSTAB  method  are  given  in  Fig.  28. 
The  main  difference  of  FLEXSTAB 's  method  compared  to  the  Koerner ' s 
technique  is  that  the  treatment  on  the  "body"  is  different. 
FLEXSTAB' s method  treats  the  body  by  a number  of  discrete  vortices 
distributed  ovi  the  mean  body  surface.  The  main  body  surface  is 
approximated  by  a number  of  rectangular  strips  with  free  stream- 
wise  surfaces  (Ref.  5).  In  addition,  FLEXSTAB's  method  was  more 
complicated  and  allows  one  to  take  into  account  the  thickness, 
the  camber  and  the  twist  effects  of  a wing  by  a distribution  of 
additional  source-sink  and  vortex  singularities. 


Figure  27.  Computational  wing-body  configuration  (Red.  f 8 J ) . 
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(e)  FLEXSTA8  PANELING  SCHEME.  (R»f.  8) 
Figure  28.  Paneling  schemes. 


2.1  Effect  of  Angles  of  Attack  on  Chordwise  Pressure 

Distribution. 

The  calculated  chordwise  pressure  distribution  on  the 
wing  is  shown  in  Fig.  29,  in  which  the  angle  of  incidence  of  the 
body  was  taken  two  degrees,  and  the  free-stream  Mach  number  was 
0.85.  The  computation  was  made  at  three  different  spanwise 
points.  The  agreement  with  FLEXSTAB's  result  was  surprisingly 
good.  The  agreement  with  the  experimental  data  was  also  satis- 
factorily good,  except  in  the  region  very  close  to  the  wing  lead- 
ing edge,  where  a sudden  pressure  drop  can  be  seen  because  of  the 
local  flow  separation.  The  case  of  the  angle  of  incidence  of 
eight  degrees  is  shown  in  Fig.  30.  The  agreement  with  FLEXSTAB's 
method  was  very  good  in  this  case  also.  The  agreement  with  the 
experimental  data  was  still  very  good  at  the  region  very  close 
to  the  wing  root  although  the  small  fluctuation  of  ACp  was  caused 
by  the  present  numerical  calculation.  However,  the  agreement 
with  data  was  not  so  good  from  the  leading  edge  to  the  mid-chord. 
The  reason  for  this  disagreement  with  data  may  be  attributed  to 
the  flow  separation  from  the  sharp  leading  edge  because  of  a con- 
siderably higher  angle  of  incidence.  A further  study  on  such  a 
situation  is  needed. 
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2.2  Flap  Effect  on  Chordwise  Pressure  Distribution. 

A further  comparison  with  the  FLEXSTAB  method  was 
made  on  the  effect  of  the  flap  of  the  wing.  The  calculated 
chordwise  pressure  distributions,  based  on  both  methods,  are 
given  in  Figs.  31  to  33.  In  these  cases,  it  is  assumed  that  the 
wing  has  zero  angle  of  attack  and  the  flap  angle  is  8.3°.  The 
following  two  examples  (Figs.  31  and  32)  show  the  chordwise  Cp 
distribution  on  the  upper  and  the  lower  surfaces  of  a wing  and 
deflected  flap  at  = 0.4  and  0.85,  respectively.  The  calcula- 
ted position  is  near  the  mid  span.  ACp/2  was  taken  for  Cp  in 
the  present  calculation  (see  Eq . 16).  Fig.  31  shows  the 

chordwise  C distribution  at  = 0.4.  No  correction  about  the 
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round  leading  edge  of  a flat  wing  was  made  in  the  present  analysis. 
Therefore,  the  Cp  close  to  the  leading  edge  of  the  wing  has 
resulted  in  a rather  poor  agreement  with  those  of  the  FLEXSTAB  and 
experimental  data.  Except  for  this  region  of  the  wing,  the  agree- 
ment of  the  Cp  distribution  calculated  by  the  present  analysis 
compared  favorably  with  data.  The  poor  agreement  with  data  near 
the  flap  hinge  line  is  attributed  to  the  local  flow  separation  on 
the  wing  because  of  a considerably  higher  flap  angle. 

Fig.  32  shows  the  chordwise  Cp  distribution  at  = 0.85. 

In  this  case,  the  values  of  Cp  from  Ref.  (8)  were  replotted  for 
eliminating  the  influence  of  the  round  leading  edge  of  the  wing 
at  6p  = 0°  as  shown  in  this  figure.  The  agreement  of  the  present 
calculation  with  that  of  the  FLEXSTAB  was  good,  in  general,  over 
the  whole  surface  of  the  wing,  and  with  the  experimental  data, 
except  the  regions  very  close  to  the  wing  leading  edge  and  the 
flap  hinge  line.  The  poor  agreement  with  data  near  the  flap  hinge 
line  is  a result  from  the  same  reason  as  mentioned  before.  The 
variation  of  Cp  with  the  free  stream  Mach  number  was  very  small 
in  this  type  of  body  configuration. 

The  calculated  chordwise  pressure  (AC  ) distribution  at 
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three  different  spanwise  positions  at  = 0.85  is  shown  in  Fig.  33. 
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Figure  31.  Surface  pressure  distribution  of  flat  wing 
(trailing-edge  6 - 8.3*,  a - 0,  M = 0.4) 
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The  agreement  on  both  methods  is  good,  however,  the  inviscid 
flow  theory  predicts  a larger  pressure  difference  than  data 
indicated  over  the  flap  surface.  Again,  such  a disagreement  is 
probably  attributed  to  the  flow  separation  on  the  upper  wing 
because  of  a considerably  higher  flap  angle  given  in  this 
example . 


2.3  Spanwise  Load  Distribution. 

The  spanwise  load  (directly  related  to  the  lift) 
distribution  was  computed  for  two  different  angles  of  attack 
cases.  And  then  they  were  compared  with  the  results  by  the 
FLEXSTAB  method.  They  are  shown  in  Fig.  34  (for  the  chordwise 
pressure  distributions,  see  Figs.  29  and  30).  It  was  assumed 
that  there  is  little  difference  between  the  lift  and  the  normal 
force  in  this  case,  or,  it  is  assumed  negligibly  small.  In  the 
case  of  the  angles  of  attack  of  two  degrees,  the  agreement  of 
the  present  calculation  with  both  the  FLEXSTAB  and  data  is 
satisfactorily  good;  this  is  apparent  from  Figs.  29a  to  c.  On 
the  other  hand,  the  poor  agreement  with  the  experimental  data 
occurred  in  the  region  from  the  mid  span  to  the  wing  tip  at 
the  angles  of  attack  of  eight  degrees.  This  can  be  also  under- 
stood by  looking  back  at  Figs.  30b  and  c.  The  satisfactorily 
good  agreement  can  be  seen  at  the  region  below  the  mid  span. 

Figs.  35  and  36  show  cne  spanwise  load  distribution 
in  the  case  of  a deflected  flap  at  Mach  number  0.4  and  0.85, 
respectively.  Nonetheless,  even  with  the  large  flap  deflections, 
both  the  theoretical  calculations  were  in  good  agreement  with 
the  experimental  data.  This  means  that  the  flap  deflection  does 
not  influence  the  spanwise  load  distribution  so  much  in  this  case 
as  the  angle  of  attack  case  did  (see  the  case  at  a = 8°  in  Fig. 
34),  because  the  part  of  the  flap  is  small  as  compared  with  the 
total  wing  surface  area.  Little  difference  in  the  compressibility 
influence  can  be  seen  in  this  wing-body  configuration  (also,  see 
Fig.  37  for  the  compressible  effect  to  the  normal  force  slope). 
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2.4  Normal  Force  Coefficient  of  Wing-Body  Combination. 

Finally,  the  computed  result  on  the  normal  force 
coefficient  has  been  compared  as  shown  in  Fig.  37a.  The  agree- 
ment: with  the  experiments  is  very  good.  The  compressibility 
effect  based  on  the  Goethert's  similarity  rule  is  again  proven 
to  be  very  satisfactory.  A further  study  based  on  Woodward's 
method  (Ref.  6)  to  extend  the  computation  into  supersonic  flow 
region  has  been  undertaken  currently.  The  computed  result  can 
be  seen  in  Fig.  37b  and  this  will  be  discussed  in  detail  in  our 
later  publication  (Ref.  13). 


* I 

_2  0.02 


e 

U 


THEORY  (PRESENT  CALCULATION) 

THEORY  (Ref.  8) 

O EXPER'  NVlRcf.  8) 


0.01 


0.2  0A  0.6  0.0 

MACH  NUMBER  (M„) 

Figure  37a.  Normal  force  coefficient. 
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SECTION  VI.  COMPARISONS  OF  POTENTIAL  FLOW  THEORIES  WITH 
DATA  ON  CRUCIFORM  SLENDER  CLIPPED  DELTA  FIN- 
BODY  COMBINATION  WITH  ELEVON  DEFLECTION. 

The  computational  model  employed  for  a cruciform  slender 
clipped  delta  fin-body  combination  is  given  in  Fig.  38.  Each 
fin  has  a small  control  fin  called  "eleven,"  and  it  is  assumed 
that  its  hinge  line  is  located  near  the  trailing  edge  of  the 
main  fin  as  shown.  Paneling  scheme  is  shown  in  Fig.  39  for  a 
deflected  elevon. 

1.  Results  and  Discussions  on  Various  Effects  of  Elevon 
Deflection  on  Fin-Body  Combination. 

.1.1  Zero  Elevon  Deflection  (a^  = 8.5°,  = 0). 

The  first  calculation  was  made  on  the  chordwise  pressure 
distribution  at  three  different  spanwise  positions.  The  angle 
of  attack  of  the  body  is  8.5  degrees  with  zero  elevon  deflection, 
and  the  computed  results  are  shown  in  Fig.  40a.  The  investigated 
Mach  number  was  zero  and  0.8,  respectively.  The  variation  of  the 
chordwise  pressure  distribution  at  various  different  spanwise 
positions  is  given.  It  is  found  that  the  influence  of  the  free 
stream  Mach  number  is  very  small.  Fig.  40b  shows  the  longitudinal 
pressure  distribution  at  three  different  azimuthal  positions  on 
the  body  corresponding  to  the  fin  in  Fig.  40a.  It  can  be  seen 
that  the  trend  of  pressure  distribution  is  very  close  to  the  one 
in  Fig.  19b,  i.e.,  the  influence  of  fin  on  the  body  is  strong  in 
the  fin  location  part  on  the  body,  and  little  influence  from  it 
outside  of  this  region. 
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Figure  40a.  Chordwise  pressure  distribution  at  three  different  spanwise 
positions  on  fin  with  zero  deflected  elevon  of  slender  fin-body 
combination  (at,  = 8.5°,  6 = 0°,  M =0  and  0.8). 
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1.2  Deflected  Elevon  (a^  =0,  6g  = 2°). 

The  effect  of  elevon  deflection  has  been  investigated 
next.  The  chordwise  pressure  distribution  on  the  fin  is  shown 
in  Fig.  41a.  The  angle  of  deflection  of  flap  is  two  degrees, 
with  the  trailing  edges  of  the  elevons  on  both  fins  down 
symmetrically,  with  no  angles  of  attack,  and  the  free  stream 
Mach  number  is  zero,  respectively.  The  influence  of  elevon  on 
the  main  fin  was  very  small  in  such  small  elevon  deflection 
angles.  However,  such  an  influence  will  be  rather  strong  as  it 
is  closer  to  the  fin  tip.  The  corresponding  longitudinal  pressure 
distribution  on  the  body  is  shown  in  Fig.  41b.  The  influence  of 
the  elevon  deflection  on  the  body  was  very  small  in  this  case. 
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Figure  41a.  Chordwise  pressure  distribution  on  fin  with  deflected 
elevon  (6e  = -2°,  <5u  = 0,  M =0). 


1.3  Combined  Effects  (a^  = 8.5°,  6g  = 2°). 

Fig.  42a  shows  the  chordwise  pressure  distribution 
on  a fin  at  an  angle  of  attack  of  the  body  of  8.5°  as  well  as 
the  symmetrical  elevon  deflection  of  2°  or  -2°,  where  the  results 
of  Figs.  40a  and  41a  were  combined  linearly.  The  corresponding 
longitudinal  pressure  distribution  on  the  body  is  shown  in  Fig. 
42b.  Thrs  result  was  made  also  by  the  linear  combination  of 
the  ones  of  Figs.  40b  and  41b. 
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gure  42b.  Longitudinal  body  surface  pressure  distribution  with  and  without  elevon 
deflection  (M  = 0,a,  = 8.5°). 


1.4  Asymmetrically  Deflected  Elevon. 

The  roll  effect  of  the  elevon  was  investigated  on 
the  cruciform  fin-body  combination  at  = 0,  and  is  shown  in 
Fig.  43.  The  angle  of  deflection  of  a pair  of  horizontal 
elevons  was  taken  as  2°  asymmetrically  in  the  negative  direction 
of  roll  in  sign  convention.  Fig.  43a(a)  shows  the  load  distri- 
bution on  the  horizontal  fin,  and  Fig.  43a (b)  shows  the  one  on 
the  vertical  fin  (i.e.,  no  deflected  elevon  case),  respectively. 
The  influence  of  the  elevon  deflection  on  the  vertical  fin  was 
not  so  large,  but  some  amount  of  induced  pressure  difference 
(it  acts  as  positive  roll;  see  Ref.  (2)  for  more  details  about 
this)  can  be  seen  from  Fig.  43a(b).  Fig.  43b  shows  the  corres- 
ponding longitudinal  Cp  distribution  on  the  body.  The  influence 
of  the  elevon  deflection  on  the  body  was  very  small  in  this  case. 
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1.5  Fin-Body  Combination  with  Angles  of  Attack  and  Yaw. 

The  pressure  distribution  on  the  body  is  shown  in 
Fig.  44  for  the  fin-body  combination  at  angles  of  incidence  of 
7.4°.  This  result  was  obtained  from  the  linear  summation  of  the 
ones  of  angles  of  attack  of  6.8°  and  yaw  angle  of  2.8°  based  on 
the  linearity  of  the  basic  equations  (see  Ref.  (2)). 
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Figure  44.  Longitudinal  pressure  distribution  on  body  with  angle 
of  attack  and  yaw  (0°elevon  deflection  angle). 


2.  Comparison  with  Experimental  Data. 

The  angle  of  incidence  of  a body  is  taken  at  8.54° 
with  an  elevon  deflection  angle  of  ±30°,  with  which  the  body 
generates  a counterclockwise  roll  facing  to  the  upstream  direc- 
tion. The  free  stream  Mach  number  is  1.62.  With  a proper 
sweptback,  as  in  this  example  calculation,  the  leading  edge  of 
the  fin  is  submerged  in  a subsonic  stream.  Because  of  the 
slenderness  of  the  fin  and  its  subsonic  leading  edge,  one 
expects  that  a better  agreement  may  be  resulted  by  employing  a 
slender  body  theory  than  the  singularity  distribution  method. 

After  a comparison  with  data,  it  indicates  indeed  that  this  is 
true.  (The  analysis  at  supersonic  speeds  by  using  "Singularity 
Method”  similar  to  that  at  subsonic  speeds  is  currently  under 
investigation.  The  method  employed  is  mainly  as  given  in  Ref.  (6).) 

Figs.  45a  to  f show  the  chordwise  pressure  distribu- 
tion computed  by  singularity  distribution  and  slender  body 
theories,  at  different  spanwise  positions  together  with  experi- 
mental measurements  (MICOM) . The  agreement  of  the  singularities 
distribution  method  by  Koerner  with  the  data  was  not  bad  in 
general  over  a fin  surface  even  though  the  flow  is  supersonic. 

The  pressure  distribution  on  the  elevon  surface  could  not  be 
compared  because  of  a lack  of  the  data  on  this  part. 

As  is  expected,  the  computed  result  by  slender-body 
theory  showed  a better  agreement  with  the  data  over  a fin.  This 
is  especially  so  on  the  lower  surface  of  the  fin,  however,  it 
cannot  calculate  the  effect  of  the  elevon  deflection,  and  a 
clipped  tip  effect.  The  Koerner' s method  estimated  a rather 
large  pressure  distribution  on  the  fin  close  to  the  hinge  line 
of  the  elevon  because  of  a very  large  elevon  deflection  angle. 

The  large  pressure  increase  is  again  attributed  to  the  flow 
separation.  This  can  be  seen  more  clearly  on  the  upper  surface 
of  the  fin  near  the  hinge  line.  The  tip  effect  caused  by  the 
leading  point  P of  the  tip  Mach  cone  may  cause  the  pressure  to 
increase.  However,  no  correction  has  been  made  in  the  present 
analysis . 
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Figure  45c.  Chordwise  load  distribution  on  fin  for  same  conditions 
as  Figure  69.*  Y = 0.752  in. 
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Figure  45d.  Chordwise  pressure  distribution  on  fin  for  same 
conditions  as  Figure  69:  Y ■>  1.653  in. 


3.  Variation  of  Surface  Pressure  and  Lift  Coefficient 
with  Free  Stream  Mach  Number. 

Various  compressibility  corrections  to  the  pressure 
coefficient  of  a sample  point  on  the  fin  are  shown  in  Fig.  46. 

The  angle  of  attack  of  8.5°  was  assumed  in  this  case.  The  com- 
puted results  used  Prandtl-Glauert  and  Karman-Tsien  rules  and 
showed  an  over-correction  as  discussed  earlier.  The  total  lift 
coefficient  of  this  particular  fin-body  combination  with  zero 
elevon  angle  is  shown  in  Fig.  47.  The  maximum  body  cross-section 
area  was  taken  as  the  reference  area  for  the  lift  coefficient 
calculation.  The  contribution  of  the  body  on  the  total  lift  co- 
efficient is  very  small.  Thus,  as  far  as  the  total  lift  coeffi- 
cient is  concerned,  the  effect  due  to  the  body  may  be  neglected. 
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Figure  46.  Comparison  of  compressibility  effect  of  present  result 
with  other  similarity  rules  (a.  ■=  8.5°). 
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SECTION  VII.  CONCLUSIONS  AND  RECOMMENDATIONS. 


The  computational  scheme  for  estimating  the  pressure  distri- 
bution on  a fin-body  combination  has  been  developed  by  the  "Sing- 
ularity Method"  based  on  Ref.  (1).  The  reasonable  results  were 
obtained  over  various  kinds  of  fin-body  configurations.  The 
comparisons  with  the  experimental  data  were  generally  good. 
Goethert's  similarity  rule  predicted  well  the  compressibility 
effect  up  to  very  high  subsonic  speed.  This  is  not  only  for  a 
moderately  large  aspect  ratio  fin  but  also  to  a slender  fin-body 
combination  at  small  angles  of  attack.  The  poor  agreement  of 
pressure  distribution  with  the  data  appeared  at  a relatively  high 
angle-of-attack  (eight  degrees  in  the  present  calculation)  as 
reported  in  the  other  method.  This  seemed  to  be  in  the  flow 
separation  on  the  lifting  surface  due  to  high  angles  of  attack, 
and  thus  further  study  on  this  will  be  needed.  The  calculated 
result  by  the  slender-body  theory  (Ref.  (2))  had  a surprisingly 
good  agreement  with  t>  e experimental  data  for  the  slender  fin- 
body  combination  even  at  supersonic  speed  as  expected,  except  in 
the  influence  region  of  the  fin  tip,  i.e.,  inside  the  tip  Mach 
cone.  The  present  analysis  has  coincided  well  with  the  slender-body 
theory,  for  slender  fin-body  configuration,  and  thus  it  may  be 
said  that  the  present  analysis  includes  the  slender-body  theory 
at  subsonic  speed. 

Although  a flat  plate  fin  case  only  has  been  investigated 
here,  it  is  not  too  difficult  to  modify  the  present  analysis  so 
that  it  can  be  applied  to  include  the  cases  of  a fin  with  camber, 
twist,  and  thickness.  Furthermore,  the  restriction  of  the  straight 
leading  and  trailing  edge  is  not  so  severe.  Therefore,  the  present 
analysis  can  be  applied  to  almost  all  the  arbitrary  plan  formed 
fin-body  combinations  with  a slight  modification  to  the  program. 
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APPENDIX  I.  DETERMINING  AERODYNAMIC  MATRICES 


1.  Determining  the  positions  of  a horseshoe  vortex  and  a con- 
trol point  on  a fin 

By  using  the  coordinates  as  shown  in  Figs.  2 and  3,  the 
positions  of  a horseshoe  vortex  (xv,  yv,  zv)  a control  point 
(xD,  yD,  zD)  on  a fin  can  be  written  as: 


= *d  = <ek + ek,t+i)/2 
xv  = ak  + ck  • (3dk,m  + dk,m+l)/4  + (yV  - bk>  • 


tan  A 


Zv  = ZD  = 0 


ak  + ck  ’ (dk,m  + 3dk,m+l)/4 
+ (ck+l  ' (dk+l,m  + 3dk+l,m+l)/4 


- ck  ‘ <dk,m  + 3dk,m+l>/4  + at 

<*v  ' bk>/<bk+l  - V 


where , 


ten  A - (ck+1  • (3dk+lm  + dk+lm+1)/4 


- ck  <3dk,m+  dk,m+l)/4 


+ ak+l  ~ ak ) /(bk+l 


The  half  spanwise  width  of  a horseshoe  vortex  is, 


h = (ek,m  * ck,£>/2 


1 


2.  Velocity  components  induced  by  a horseshoe  vortex 


The  vertically  intercepted  point  S ( £ , 3n . BZV)  between  a 
bound  vortex  line  and  any  control  point  P(xD,ByD  BZD)  as  shown 
in  Fig.  48  can  be  written  as: 


C 


Bn 


xv  + xD  tan2  p + B • (yD~yv)  tan  p 

' o 

tan  p + 1 
2 

eyD  + 6yv  tan  p + (xD-xv)  tan  p 
— 

tan  p + 1 


S 


where , 

tan  p = tanX/B 


fy  using  Biot-Savart  law,  the  velocity  (Wfi)  induced  at  P by  a 
bound  vertex  £,  can  be  written  as: 


T (cos  y + cos  <f>) 
= - - — 

4tt  7 A2  + Z2 


where 

± / B2  + C2 

cos  y = — = — - - ■ - 

^ y A2  + B2  + C2  + Z2 


\ 


cos  y > 0 for  n > y^  - h 
cos  Y < 0 for  n < yv  - h 


cos  4> 


RS 

KF 


± y d2  + e2 

yA2  + D2  + E2  + Z2 


cos  ^ 

< 0 

for 

n 

> 

A 

+ 

h 

COS  <t> 

> 0 

for 

n 

A 

< 

+ 

h 

(23) 


(24) 
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(SF')2  = (xD  - O2  + B2 
3 • (n  - yv  + h) 


<n  - yr,)' 


£ - xv  + Bh  tan  y 


3 • (yv  + h - n) 


= Xy  + 0h  tan  y - £ 

= 6 • (zv  - zD) 


thus,  the  velocity  components  (ug,  Vg,  wg)  of  Wg  can  be  written 


Un  = - Wt>  sin  e cos  y 


VB  = WB  s*n  e s^-n  y 


wg  = - Wg  cos  e 


where , 


sin  e 


2 2 
+ Z 


cos  e = 


A2  + Z2 


cos  e > 0 for  xD  > £ \ 

cos  c < 0 for  Xg  < ^ 

The  velocity  (Wp)  induced  at  P due  to  a port  free  vortex  as  shown 

in  Fig.  49  can  be  written  as: 


P Att 


1 + cos  k 


G2  + Z2 
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where 


cos  k 


H 


+ Zi 


G = e • (yD  - yv  + h) 

H = xd  _ + tan  P 


c 


thus,  the  velocity  components  (u  v w ) of  W can  be  written 
as : f F 


up 

0 

VP  = 

-WP 

sin  v 

li 

3°* 

- WP 

cos  V 

where , 

N 


sin  v = 


p'x.2 


G + Z‘ 
G 


COS  V = 


/ G2  + Z‘ 


(2 


Similarly,  the  velocity  (Wg)  induced  at  P due  to  a starboard  free 
vortex  as  shown  in  Fig.  50  can  be  written  as: 


WS  = - 


1 + cos  a 


Air  / 12  + Z2 


\ 


where, 

cos  a = 


/ 12  + J2  + Z2 


I " B (yv  + h - yD) 

" - xv  - Sh  tan  )i 


(21 


J 


1 


COS  T 


/r*”+  z2 


(29) 


Therefore,  the  induced  velocity  components  (u,  v,  w)  due  to  a 
complete  horseshoe  vortex  is  considered  as  a linear  summation 
of  the  three  vortices  mentioned  above,  and  can  be  written  as: 


u 


UB  + Up  + Ug 


v = 


VB  + Vp  + Vg 


w = 


Wg  + Wp  + Vg 


OO) 


or, 


T cos  y + cos  ip  „ j* « 

i n — - Z cos  y = PT 

Att  Az  + Zz 


■ x f 

4 IT  l 


cos  Y + cos 


A2  + Z2 


1 + cos  k 
sxn  y - — ^ y~ 


+ V 


COS  0+1 


— 5 ’ 

Iz  + Z‘ 


J Z = Qf 


w 


r r cos  y + cos  $ A 1 + cos  k 

£ l ' ~ '+“?  ' ' G^  + Z2 

cos  0 + 1 T j 


I2  + zl 


Rr 


The  first  term  in  w;  " £or  XD  > £ 

+ for  xc  < f, 


(31) 
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3.  Aerodynamic  Matrices. 

For  a planar  fin-body  configuration  geometry,  the  effect 
of  the  horseshoe  vortex  on  the  opposite  side's  fin  and  the  image 
vortex  can  be  treated  in  the  same  way  as  that  mentioned  above. 
Thus,  total  velocity  components  at  any  control  point  induced  by 
four  horseshoe  vortices  (i.e.,  two  on  a pair  of  .tins,  and  two 
inside  a body)  as  shown  in  Fig.  51)  can  be  written  as: 


r /p  (SS)  p (PS)  p.(SS)  + p i (PS)  v 
ul  - l (Pij  + Pij  + pij  + Fij  ■ ) rj 

j-1 


l «USS> 


Q • + Q ! V 


.(PS) 


l (RijSS)  + + RijSS)  - R:fs)  ) r . 


due  to  a pair 
of  fins  1 vortices 


due  t(>  image  vortices 
inside  a body 


where,  the  superscript  (SS)  indicates  that  as  well  vortex  as 
control  point  are  situated  on  the  starboard  side  (see  Ref.  (7)  ), 
and, 

Pf(PS)  _ p(SS)  etc>>  with  yD  = -yD  , on  the  opposite  fin 

2 

P:(PS)  = Pijss)  etc.,  with  y = y-  for  0 < y < a,  inside  a body 

It  has  been  assumed,  for  simplicity,  that  the  bound  vortex  of 
an  image  has  the  straight  line  instead  of  the  actual  curved  line. 

By  the  relation  of, 

* = EXT  * (33) 
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APPENDIX  2.  COMPUTER  PROGRAM  FOR 
„ DETERMINING  PRESSURE 

DISTRIBUTION  ON  FIN-BODY 
COMBINATION 

2.1  Program  Description 

A computer  program  has  been  developed  to  calculate  the 
pressure  distribution  and  aerodynamic  characteristics  of  fin  and 
fin-body  combination  in  subsonic  flow.  The  program  is  written 
in  FORTRAN  IV,  and  designed  for  IBM  360  computer,  however,  it  is 
easily  adaptable  to  other  computers  with  minor  modifications. 


2.2  Program  Structure 

The  computer  program  consists  of  one  main  overlay  program 
and  thirty -one  subroutines. 

Main  Program 

The  brief  outline  of  ti  j Main  program  is  shown  in  Fig.  52. 

The  iteration  scheme  is  used  only  for  a fin-body  combination  to 
solve  a fin-body  interaction  result. 

A considerably  detailed  expression  of  the  Main  program  is 
shown  in  Fig,  53.  The  Main  program  controls  almost  all  subroutines. 
The  effect  of  the  horizontal  fin  on  the  vertical  fin  is  computed 
in  the  rear  part  of  this  Main  program. 


Figure  5?..  Outline  of  Main  Program. 


Aerodynamic*  Performance  of  Planar  Fin  or 
Fin-3ody  Combination  Must  Be  Computed  at  First. 


Figure  53.  Flow  Diagram  for  MAIN  PROGRAM. 
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Call  BLLT 
VBIA 


Call  ARDC 


'JRAD' 


Call  PVMX 
ARQQ 
ARQV 
ARBF 


IRPT  = 1 


JRAD 


Call  BCF  (BC,  NO) 
VDF  (GAMM) 
PDB  (XCP) 

SLD 

XLF  (CL) 

PMS  (CMS) 


Figure  53.  (Continued) 


Computation  for  Aerodynamic  Performance  of 
Planar  Fin  or  Fin-Body  Combination  Terminates. 
Subsequently,  Effects  of  Horizontal  Fin  on 
Vertical  Fin  is  Computed. 


Call  PVMX 
ARBF 

DINV  (AR,  BC,  NO) 
PDB  (XCP) 

SLD 

XLF  (CL) 

Write  XCP 
CLS 
CLSW 
CLSF 
CL 


Write  CPB 


Figure  53.  (Concluded) 


r 
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Subprogram  SUBROUTINE  ARDC 

In  this  subprogram,  the  aerodynamic  matrices  of  a fin  to  fin 
will  be  computed.  The  subprograms  PVB,  PVP,  and  PVS  are  called 
here. 


In  this  subprogram, 
NN  = 2,  no  body 


where , 


■ 4,  with  body  (see  Fig.  54) 

AR  (I,J)  - SAR  + ASYN  * PAR  + SIAR  + ASYN  * PIAR 
influenced  panel 
influencing  panel 

for  NN  = 1 

- 2 
« 3 
- 4 


I. 

J. 
SAR. 
PAR. 

SIAR. 

PIAR. 


WB  + WP  + WS 


(see  Fig.  54) 


ASYN. 


- ASYW 

- ASYF 


for  main  fin  or  wing 
for  flap  or  elevon 


Exchange  of  the  Cards  (except  Data  cards) 

In  the  list  of  the  subroutines  BGMY  and  GVS1,  the  cosine 
paneling  is  shown.  A few  cards  must  be  exchanged  by  the  suitable 
ones  when  a user  wants  the  equidistant  paneling.  They  are  shown 
in  Table  1.  A clear  example  of  the  spanwisely  equidistant,  chord- 
wisely  cosine  paneling  on  a fin,  and  longitudinally  cosine  panel- 
ing on  a body  is  shown  in  Figs.  28(a)  and  (b) . 
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TABLE  1.  EXCHANGE  OF  THE  CARDS 


2.3  List  of  Subprograms 

All  the  Subprogram  SUBROUTINES  and  their  purposes  are 
shown  in  Table  2.  Most  of  them  are  controlled  by  the  Main  program. 

2.4  List  of  Symbols 

The  symbols  and  their  meanings  used  in  the  Main  program 
and  SUBROUTINES  are  shown  in  Table  3 . 


TABLE  2 

LISTS  OF  SUBPROGRAMS 


NAME 


PURPOSE 


Mi 


NP 


GVS  2 (SX,  SY,  TMU) 


PVB  (UB,  VB,  WB) 
PVP  (UP,  VP,  WP) 


PVS  (US,  VS,  WS) 
BCF  (XB,  N3) 


■lRDC 


Defining  location  of  control  points  on 
the  fin,  and  half  width  of  local  bound 
vortex. 


Defining  body  control  points  (which  also 
equals  to  location  of  body  sources) . 

Naming  location  of  total  control  points 
on  fin-body  combination. 


Computing  location  and  slope  of  body 
image  vortex. 

Computing  interception  at  S.  (See  Fig. 
48.) 


Computing  normalized  perturbation  velo- 
city components  induced  by  bound  vortex. 

Computing  normalized  perturbation  velo- 
city component  induced  by  port  free  vor- 
tex. 


Computing  perturbation  velocity  component 
induced  by  starboard  free  vortex. 

Defining  initial  angle  of  attack  of  fin, 
and  induced  angle  of  attack  due  to  angle 
of  incidence  of  body  (when  it  is  not  zero) 
and  distributed  sources  on  body. 

Computing  aerodynamic  matrices  of  fin 
only. 
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Table  2.  Li/Sts  of  SUBPROGRAMS  (Continued) 


NAME 

PURPOSE 

VDF  (GAMM) 

Compute  circulation  strength  (i.e., 
solve  tho  simultaneous  linear  eqs . ) • 

GVS  3 

Computing  local  fin  panel  area,  refer- 
ence fin  area,  and  reference  chord. 

XLC 

Computing  local  chord  of  fin. 

BLLT 

Computing  local  chord  length  of  fin 
provided  that  fin  was  extended  into 
body. 

DINV  (D,  E,  N) 

Solving  the  simultaneous  linear  equation, 

PDB  (XCP'> 

Computing  pressure  distribution  (ACp)  on 
fin. 

SLD 

Computing  spanwise  lift  distribution 
based  on  local  fin  chord. 

XLF  (CL) 

Computing  lift  coefficient  of  fin  based 
on  reference  area  (RFAREA) . 

PMS  (CMS) 

Computing  spanwise  moment  distribution 
around  257.  local  chord  based  on  local 
chord  CLL  and  XAX. 

PVMX 

Computing  perturbation  velocity  compo- 
nent on  body  due  to  unit  circulation  of 
fin. 

VLNB 

Computing  perturbation  velocity  compo- 
nents on  body. 

BPD 

Computing  pressure  distribution  (Cp)  on 
body. 

ARQQ 

Computing  pressure  distribution  on  body 
due  to  unit  source  strength. 
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Table  2.  Lists  of  SUBPROGRAMS  (Continued) 


NAME 

PURPOSE 

BLD  (BCL) 

Computing  pressure  distribution  (ACp)  on 
body. 

BSS 

Computing  source  strength  on  body. 

ARQV 

Computing  body  source  strength  for  m^ 
iteration. 

WBS  (XW,  Nl,  N2) 

Computing  normalized  downwash  on  fin  in- 
duced by  body  source. 

ARBF 

Computing  normalized  downwash  on  fin  in- 
duced by  body  unit  source. 

BSLD 

Computing  spanwise  lift  distribution  on 
body  based  on  local  chord  of  fin,  sup- 
posed to  be  extended  into  body. 

VBIA 

Computing  angle  of  attack  on  fin  induced 
by  angle  of  incidence  of  body. 

XLIFT 

Computing  lift  coefficient  of  total  fin- 
body  combination. 

TABLE  3 


LIST  OF  SYMBOLS 


Control  parameter.  Integer  inside  of  ( ) denotes  dimensions 


SYMBOLS 


MEANINGS 


BD  (56) 


XW  (56) 


GAMM  (55) 
XREF  (10) 
XCP  (55) 


CMS  (10) 


BCL  (36) 


CPB  (72) 

AR  (56,  56) 


AQ  (72,  72) 


ASV  (72,72) 


(Equivalent  to  XB  in  SUBPROGRAM  BCF. ) 
Apparent  angle  of  attack  of  local  panel. 

Normalized  component  of  induced  velocity 
on  local  panel  due  to  body  source. 

Circulation  strength  of  local  panel 

X-coordinate  for  257.  local  chord  of  fin. 

Pressure  difference  (AC  ) on  local  panel 
of  fin. 

Spanv/ise  moment  distribution  around  257. 
local  chord,  based  on  local  chord  of 
fin. 

Pressure  difference  (ACp)  on  body. 

Radius  of  body  (=  constant  in  present 
analysis) . 

Pressure  coefficient  (Cp)  on  body. 

Aerodynamic  influence  matrix  for  fin 
only. 

Perturbation  velocity  component  (x-direc 

tion)  on  body  due  to  body  source  ((Cp)  = 

- 2 (u  + u ) /B2 , see  SUBPROGRAM  BPD) . B 
s 

1 £ l-cos (0^ -9 v)  ] BaAS 


2n  [(xp-xv)2  + B2  { 1-cos (0^ -0 v)  } a2] 
(See  Eq.  (10).) 


Table  3.  'List  of  Symbols  (Continued) 


SYMBOLS 


MEANINGS 


AFN 

(72,55) 

MB 

NB 

NUB 

PSA 

(55) 

BREF 

(5) 

XDV 

(55) 

YVI 

(55) 

WHI 

(55) 

TLMDI  (55) 

NBF, 

NBM,  NBB 

Coefficient  for  determining  induced 
velocity  (w-component)  on  fin  due  to  unit 
source  strength  on  body. 

Circumferential  division  number  along 
half  body. 

= NBF  + NBM  + NBB . 

= NB  x MB;  total  panel  numbers  on  half 
body. 

Local  panel  area  (includes  compressible 
effect) . 

Local  chord  of  fin,  provided  extending 
into  body. 

“ XG  - SV;  1/2  local  panel  chord  through 
its  center  in  present  analysis. 

Y - location  of  center  of  bound  vortex 
for  body  image. 

Half  width  of  body  image  vortex. 

Slope  (or  gradient)  of  bound  vortex  of 
body  image. 

Longitudinal  division  numbers  on  body  for 
forward  (ahead  of  L.  E.  of  fin  location), 
middle  (at  fin  location) , and  backward 
(aft  of  it),  respectively. 

See  Fig.  56(a). 

Angle  of  attack  of  body  (degrees) . 


Table  3.  * List  of  Symbols  (Continued) 


SYMBOLS 


MEANINGS 


XBW  (55) 


AREAB  (96) 
UBDY  (72) 


VN  (72) 


QL  (72) 


Additive  angle  of  attack  of  local  panel 
of  fin  due  to  body  angle  of  attack. 

Control  integer.  See  Figs.  52,  53. 

Local  panel  area  on  body  surface. 

Induced  velocity  (u-component)  at  control 
point  on  body. 

Induced  velocity  component  normal  to  body 
surface  due  to  circulation  of  fin. 

Source  strength  per  unit  body  surface 
area . 

Control  variables.  See  2.2  Input  Data  Cards 


CLL  (55) 


AFLAP 


Local  chord  length  ( = XC3  - XC1) . 

Load  distribution  on  body. 

*=  MB/2. 

Chordwise  division  number  for  flap  or 
elevon. 

*=  N + NFLP. 

«=  M x N. 

Flap  or  elevon  deflection  (degrees) . 
Control  integer.  (See  Figs.  52,  53). 
See  Fig.  55(a) . 

Spanwise  division  number  for  main  fin. 
Chordwise  division  number  for  main  fin. 


Table  3.  - List  of  Symbols  (Continued) 


SYMBOLS 


ALFAF 


XCl 

(55) 

XC2 

(5) 

XC3 

(5) 

ABU 

(72, 

55) 

ABV 

(72, 

55) 

ABW 

(72, 

55) 

XL 

(10) 

XG 

(55) 

YG 

(55) 

ZG 

(55) 

XS 

(72) 

TS 

(127) 

Free  stream  Mach  number. 

Angle  of  attack  for  fin  (degrees) . 

M x NT  (total  panel  number  for  half  fin) 
Constant  (=  3.14159). 

- h - 

(-  v 2) 

Location  in  x-direction  of  intercept 
points.  See  Fig.  55(c). 


Induced  velocity  components  at  control 
point  on  fin,  due  to  circulation  of  fin. 


XC2  - XCl.  See  Fig.  55(c). 

“ NO  + NUB,  total  panel  number. 
Location  of  control  point  on  fin. 


Location  of  source  distribution  on  body 
(=  location  of  control  point  on  body) . 

Angular  position  for  source  (thus,  control 
point)  on  body  (radian) . I 


Table  3.  'List  of  Symbols  (Continued) 


SYMBOLS 


XD  (127) 

YD  (127) 

ZD  (127) 

XV  (55),  Y V 


MEANINGS 


Location  of  control  point  of  local  panel 
l of  fin-body  combination. 


WH  (55) 
TLMD  (55) 


YF 


RFAREA 

XAX 


(55) , ZV  (55)  Location  of  middle  of  bound  vortex  on 
local  panel  of  fin. 

Half  width  of  local  panel  on  fin. 

Slope  of  local  bound  vortex.  See  Fig. 
55(b). 

Y - location  of  local  panel  edges.  See 
Fig.  55(b). 

Reference  area  for  C, 


DT 


*IRC 
DX  (96) 
CLS  (5) 


CLSW  (5) 
CLSF  (5) 
CLBS  (5) 


'L- 


Maximum  chord  length  of  fin,  reference 
chord. 


Circumferential  length  of  local  body 
panel . 

Control  integer.  See  Fig.  53. 
Longitudinal  length  of  local  body  panel. 
Spanwise  lift  distribution  for  wing  + 


fin  based  on  CLL. 


Spanwise  lift  distribution  for  wing, 
based  on  CLL. 


Spanwise  lift  distribution  for  fin, 
based  on  CLL. 


Spanwise  lift  distribution  on  body,  based 
on  BREF. 


i 
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Table  3. 

t 

SYMBOLS 

CLWING 

CLFIN 

CLBODY 

TOTCL 


List  of  Symbols  (Continued) 


MEANINGS 

CL  for  wing  or  main  fin,  based  on  refer- 
ence area  (RFAREA) 

CL  for  flap  or  elevon,  based  on  reference 
area  (RFAREA) . 

CL  for  body,  based  on  RFAREA. 

Total  CL  (-  <CL)utng  + <CL)fin  + <CL)body  ) 


; 


i 


2.5  INPUT  DATA  Cards  " 

The  fin-*body  geometry  must  be  given  by  the  1st  and  2nd  cards. 


Column 

Symbols 

Explanation 

1st 

Card 

FORMAT  (8F10.5) 

1-10 

XI 

Coordinate  of  the  corner 

11  - 20 

Y1 

points  of  main  fin  or  wing. 

21  - 30 

X2 

(See  Fig.  55(a).)  If  no 

31  - 40 

Y2 

flap  nor  elevon,  X3 , Y3,  X4, 

41  - 50 

X3 

(or  X5) 

Y4  must  be  the  same  values 

51  - 60 

Y3 

(or  Y5) 

as  X5,  Y5,  X6,  Y6,  respec- 

61 - 70 

X4 

(or  X6) 

tively. 

71  - 80 

Y4 

(or  Y6) 

2nd 

Card 

FORMAT  (4F10.5) 

1-10 

X5 

Coordinate  of  the  corner 

11  - 20 

Y5 

points  of  a flap  or  an  ele- 

21 - 30 

X6 

von. 

31  - 40 

Y6 

Panel 

number  of  a fin 

must  be  given  by  the  3rd  card. 

3rd 

Card 

FORMAT  (3110) 

1-10 

M 

Spanwise  panel  number  of 
main  fin.  (See  Fig.  55(b).) 

11  - 20 

N 

Chordwise  panel  number  of 
main  fin. 

21  - 30 

NFLP 

Chordwise  panel  number  of 

flap  or  elevon. 

M must  be  less  than  or  equal  to  five.  The  sum  of  N and  NFLP 
must  be  less  than  or  equal  to  eleven. 

Panel  number  of  a body  must  be  given  by  the  4th  card. 

4th  Card  FORMAT  (4110,  2F10.5) 


1 - 

10 

MB 

Circumferential  panel  number. 
(See  Fig.  56(a).) 

11  - 

20 

NBF 

i 

! Longitudinal  ( for  fronL 

part 

21  - 

30 

NBM  | 

panel  number  | icn 

part 

31  - 

40 

NBB 

1 v for  aft 

part 
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4th  Card  (Continued) 


FDL 

BDL 


Front  part  length/ (X3-X1) . 
Aft  part  length/ (X3-X1) . 


41»  - 50 
51  - 60 

MB  must  be  six.  Total  sum  of  NBF,  NBM,  and  NBB  must  be  less 
than  or  equal  to  twelve. 


Aerodynamic  Input  must  be  given  by  the  5th  card. 
5th  Card  FORMAT  (6F10.5) 


1 - 

10 

XM 

Free  stream  Mach  number. 
(See  Fig.  57.) 

11  - 

20 

ALFAF 

Angle-of-attack  of  main  fin 
or  wing. 

21  - 

30 

ALFAB 

Body  Angle-of-attack. 

31  - 

40 

AFLAP 

Angle-of-attack  of  flap  or 
elevon. 

41  - 

50 

ASYW 

= 1.  for  symmetrical  flow 
of  main  fin  or  wing; 

■ -1.  for  asymmetrical  flow 
of  main  fin  or  wing. 

51  - 

60 

ASYF 

= 1.  for  symmetrical  flow 
of  flap  or  elevon. 

■ -1.  for  asymmetrical  flow 
of  flap  or  elevon. 

XM  must  be 

less  than 

1.0.  i 

Only  when  ALFAB  equals  to  zero, 

ASYW  and/or  ASYF  can  take 

-1.0. 

Nonetheless  that  ASYW  and/or  ASYF 

equal  to  1.0  or 

-1.0,  the 

same  amount  of  angle-of-attack  of  a pair 

of  fins  must  be  taken,  thus  the  case  that  the  right  hand  side  of 
a pair  of  fins  has  angle-of-attack  of  5°,  and  the  left  has  3°, 
cannot  be  allowed. 
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Figure  56b.  Body  paneling  and  body  control  point. 
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2.6  Program  Output 
» 

One  example  of  the  computer  print  output  is  shown  in 
Table  4.  The  corresponding  fin  is  the  following: 

Fin-body  configuration cylindrical  body  with 

rectangular  fin 

Fin  aspect  ratio 6 

Maximum  span/body  diameter..  6 

Fin  panel  number 25 

Body  panel  number 60 

Total  fin-body  panel  number.  85 

Free  stream  Mach  number 0.4 

Angle  of  incidence  of  body. . +2(°) 

No  cant  angle  of  fin,  no  flap 

Remark:  Note  that  WING  CL,  FIN  CL,  BODY  CL,  and  TOTAL  CL 
in  Output  are  half  values  of  the  correct  ones,  thus  these  values 
have  to  be  multiplied  by  two  to  obtain  correct  CL  based  on 
reference  area  (also,  see  SUBROUTINE  XLIFT) . CL  for  fin  only  is 
also  in  the  same  thing. 
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